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Abstract. It is well known that li(a;) > tt(x) (i) up to the (very large) 
Skewes' number x\ ~ 1.40 X 10 316 pQ. But, according to a Littlewood's 
theorem, there exist infinitely many x that violate the inequality, due to 
the specific distribution of non-trivial zeros 7 of the Riemann zeta function 
C(«), encoded by the equation \\{x) - tt(x) ^^[1 + 2E 7 SMxMM] 
If Riemann hypothesis (RH) holds, (i) may be replaced by the equivalent 
statement li[i/>(x)] > iz(x) (ii) due to Robin [2]. A statement similar to (i) 
was found by Chebyshev that 7r(x;4,3) — 7r(x;4, 1) > (iii) holds for any 
x < 26861 [3] (the notation n(x; k, I) means the number of primes up to 
x and congruent to I mod k). The Chebyshev 's bias (iii) is related to the 
generalized Riemann hypothesis (GRH) and occurs with a logarithmic den- 
sity 0.9959 [3|- In this paper, we reformulate the Chebyshev's bias for a 
general modulus q as the inequality B(x;q,R) — B(x;q,N) > (iv), where 
B(x; k, I) = li[0(fc) * ip{x; k, I)] — </>(fc) * k(x\ k, I) is a counting function intro- 
duced in Robin's paper [2] and i?( resp. TV) is a quadratic residue modulo q 
(rcsp. a non-quadratic residue). We investigate numerically the case q = 4 
and a few prime moduli p. Then, we proove that (iv) is equivalent to GRH 
for the modulus q. 



1. Introduction 

In the following, we denote -k(x) the prime counting function and 7r(x; q, a) the 
number of primes not exceeding x and congruent to a mod q. The asymptotic law 
for the distribution of primes is the prime number theorem 7r(x) ~ i^^- Corre- 
spondingly, one gets [H eq. (14), p. 125] 

(1.1) ^( x;g , a )^_^ 

that is, one expects the same number of primes in each residue class a mod q, if 
(a, q) = 1. Chebyshev's bias is the observation that, contrarily to expectations, 
tt(x; q, N) > -k{x] q, R) most of the times, when N is a non-square modulo q, but R 
is. 

Let us start with the bias 

(1.2) 5{x,A) :=7r(x;4,3)-7r(ar;4,l) 

found between the number of primes in the non-quadratic residue class iV = 3 
mod 4 and the number of primes in the quadratic one R = 3 mod 4. The values 
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<5(10", 4), n < 1, form the increasing sequence 

,4091295 = {1, 2, 7, 10, 25, 147, 218, 446, 551, 5960,...}. 

The bias is found to be negative in thin zones of size 

{2, 410, 15 358, 41346, 42 233 786, 416 889 978, . . .} 

spread over the location of primes of maximum negative bias [5] 

{26861, 623 681, 12 366 589, 951 867 937, 6 345 026 833, 18 699 356 321 . . .}. 

It has been proved that there are are infinitely many sign changes in the Cheby- 
shev's bias (jl.2[) . This follows from the Littlewood's oscillation theorem [H[7] 



(1.3) S(x,A) :=fi± ( — log 3: 



x^ 2 
\osx 



A useful measure of the Chebyshev's bias is the logarithmic density [3l [6j [8] 

(1.4) d(A) = lim — !— V - 

° aeA,a<x 

for the positive A + and negative A - regions calculated as d(A + ) = 0.9959 and 
d(A-) = 0.0041. 

In essence, Chebyshev's bias 5(x, 4) is similar to the bias 

(1.5) 5(x) :=H(x) -tt(x). 

It is known that 6(x) > up to the (very large) Skewes' number x\ ~ 1.40 x 10 316 
but, according to Littlewood's theorem, there also are infinitely many sign changes 

of<y(x) 0. 

The reason why the asymmetry in (jl.5J) is so much pronounced is encoded in the 
following approximation of the bias [3j I9] 1 ! 

(1.6) «(*)~^(l + 2£^ 7 ^* + a ^ 



logx ^ ^ ^l/^ + l 2 J 

where a 7 = cot _1 (27) and 7 is the imaginary part of the non-trivial zeros of the 
Riemann zeta function £( s )- The smallest value of 7 is quite large, 71 ~ 14.134, 
and leads to a large asymmetry in (|1.5p . 

Under the assumption that the generalized Riemann hypothesis (GRH) holds 
that is, if the Dirichlet L-function with non trivial real character K4 

(1.7) L( S}K4 ) = J2 [ ) 



n>0 



(2n + l) s ' 



has all its non-trivial zeros located on the vertical axis 3?(s) = 1, then the formula 
(|1.6[) also holds for the Chebyshev's bias 5 (x, 4). The smallest non-trivial zero of 
L{s, K4) is at 7! ~ 6.02, a much smaller value than than the one corresponding to 
C(s), so that the bias is also much smaller. 



The bias may also be approached in a different way by relating it to the second order Landau- 
Ramanujan constant |10] , 
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A second factor controls the aforementionned assymmetry of a L-function of real 
non-trivial character k, it is the variance 

< L8) "W = Ei7iT7- 

7>0 ' ' 

For the function £(s) and L(s, k&) one gets V = 0.045 and V=0.155, respectively. 

Our main goal. In a groundbreaking paper, Robin reformulated the unconditional 
bias (jl.5l) as a conditional one involving the second Chebyshev function if)(x) = 

(1.9) The equality 5'(x) := li[ijj(x)] — tt(x) > is equivalent to RH. 

This statement is given as Corollary 1.2 in |12) and led the second and third author 
of the present work to derive a good prime counting function 

3 

(1.10) n(x) = Y / Kn)m^) 1/n }- 

n=l 

Here, we are interested in a similar method to regularize the Chebyshev's bias 
in a conditional way similar to (|1.9[> . In [2J, Robin introduced the function 

(1.11) B(x; q, a) = \i[<j){q)^{x\ q, a)} - (f>(q)Tr(x; q, a), 

that generalizes (|1.9[) and applies it to the residue class a mod q, with ip( x y 9: a ) the 
generalized second Chebyshev's function. Under GRH he proved that [2j Lemma 
2, p. 265] 

(1.12) B(x;q,a) = fi± 
that is 



log 2 x 



(1.13) The inequality B(x;q,a) > is equivalent to GRH. 

For the Chebyshev's bias, we now need a proposition taking into account two 
residue classes such that a — N(& non-quadratic residue) and a — R (a quadratic 



one 



Proposition 1.1. Let B(x; q, a) be the Robin B-function defined in (|1.11[) . and R 
(resp. N) be a quadratic residue modulo q (resp. a non-quadratic residue), then 
the statement S'(x,q) := B(x;q,R) — B(x;q,N) > 0, Vx (i), is equivalent to GRH 
for the modulus q. 

The present paper deals about the numerical justification of proposition 11.11 in 
Sec. [5] and its proof in Sec. [3J The calculations are performed with the soft- 
ware Magma [13] available on a 96 MB segment of the cluster at the University of 
Franche-Comte. 

2. The regularized Chebyshev's bias 

All over this section, we are interested in the prime champions of the Chebyshev's 
bias S(x,q) (as defined in (|1.2I) or (|2.3p . depending on the context). We separate 
the prime champions leading to a positive/negative bias. Thus, the n-th prime 
champion satisfies 

(2.1) 6 (e) (x n ,q) =en, e = ±1. 
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Figure 1. The normalized regularized bias 6' (x , 4) / y/x versus 
the Chebyshev's bias S(x,A) at the prime champions of 5(x,4) 
(when S(x, 4) > 0) and at the prime champions of — S(x, 4) (when 
S(x,A) < 0). The extremal prime champions in the plot are 
x = 359327 (with S = 105) and x = 951867937 (with 5 = -48). 
The curve is asymmetric around the vertical axis, a fact that re- 
flects the asymmetry of the Chebyshev's bias. As explained in the 
text, a violation of GRH would imply a negative value of the regu- 
larized bias 5' (x, 4). The small dot curve corresponds to the fit of 
5'{x,4)/^fx by 2/ log x calculated in Sec. 3. 

We also introduce a new measure of the overall bias 6(g), dedicated to our plots, as 
follows 

(2.2) % ) = £^^- 

Xn 

Indeed, smaller is the bias lower is the value of b(q). Anticipating over the results 
presented below, Table [1] summarize the calculations. 

Table 1. The new bias (|2.2j) (column 2) and the standard loga- 
rithmic density (|1.4j) (column 3). 



modulus q 


bias b(q) 


log density d(A + ) 


first zero 71 


4 


0.7926 


0.9959 [3] 


14.134 


11 


0.1841 


0.9167 [3] 


0.2029 


13 


0.2803 


0.9443 [3] 


3.119 


163 


0.0809 


0.55 [9] 


2.477 
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Chebyshev's bias for the modulus q = 4. As explained in the introduction, 
our goal in this paper is to reexpress a standard Chebyshev's bias 5(x,q) into 
a regularized one S'(x, q), that is always positive under the condition that GRH 
holds. Indeed we do not discover any numerical violation of GRH and we always 
obtains a positive i5'(ie,(7). The asymmetry of Chebyshev's bias arises in the plot 
S vs 5' , where the fall of the normalized bias is faster for negative values of 5 
than for positive ones. Fig. [1] clarifies this effect for the historic modulus q = 4. 
We restricted our plot to the champions of the bias 5 and separated positive and 
negative champions. 

Chebyshev's bias for a prime modulus p. For a prime modulus p, we dehne the 
bias so as to obtain an averaging over all differences tt(x;p,N) — n(x]p, i?), where 
as above N and R denote a non-quadratic and a quadratic residue, respectively 



where ( is the Legendre symbol. Correspondingly, we define the regularized bias 
as 



Proposition 2.1. Let p be a selected prime modulus and S'(x,p) as in (|2.4| then 
the statement S'(x,p) > 0, Vx, is equivalent to GRH for the modulus p. 

As mentioned in the introduction, the Chebyshev's bias is much influenced by 
the location of the first non-trivial zero of the function L(s, K q ), n q being the real 
non-principal character modulo q. This is especially true for L(s, K163) with its 
smaller non-trivial zero at 7 ~ 0.2029 [9]. The first negative values occur at 
{15073,15077,15083,...}. 

Fig. [2] represents the Chebyshev's bias 5' for the modulus q — 163 versus the 
standard one 5 (thick dots). Tha asymmetry of the Chebyshev's bias is revealed at 
small values of \8\ where the the fit of the regularized bias by the curve 2/ log a; is 
not good (thin dots). 

For the modulus q — 13, the imaginary part of the first zero is not especially 
small, 71 ~ 3.119, but the variance (|1.8|) is quite high, V(k-i3) ~ 0.396. The first 
negative values of S(x, 13) at primes occur when {2083, 2089, 10531, . . .}. Fig. [3] 
represents the Chebyshev's bias 5' for the modulus q = 13 versus the standard one 
S (thick dots) as compared to the fit by 2/ log a: (thin dots). 

Finally, for the modulus q — 11, the imaginary part of the first zero is quite 
small, 7! ~ 0.209, and the variance is high, V(k-ii) ~ 0.507. In such a case, as 
shown in Fig. 31 the approximation of the regularized bias by 2/ log a; is good in 
the whole range of values of x. 



For approaching the proposition 11.11 we reformulate it in a simpler way as 

Proposition 3.1. One introduces the regularized couting function %'(x;q,l) := 
n(x;q,l) — ip(x;q,l)/\ogx. The statement 7t'(x; q, N) > ir'(x; q, R), \fx (ii), is equiv- 
alent to GRH for the modulus q. 



(2.3) 





(2.4) 




3. Proof of proposition II. II 
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Figure 2. The normalized regularized bias 6' (x, 163)/^/x ver- 
sus the Chebyshev's bias 5(x, 163) at all the prime champions of 
\S(x, 163)| [from \S(x, 163) | > 74 the bias is S(x, 163) < negative], 
superimposed to the curve at the prime champions of — 8(x, 163) 
(when d(x, 163) < 0). The extremal prime champions in the plot 
are x = 68491 (with S = 74) and x = 174637 (with 5 = -86). 
The asymmetry is still clearly visible in the range of small values 
of |(5| but tends to disappear in the range of high values of\6\. The 
small dot curve corresponds to the fit of 5'(x, 163) /y/x by 2/ log a; 
calculated in Sec. 3. 



Proof. First observe that proposition 11.11 follows from proposition 13.11 This is 
straightforward because according to [2] p. 260], the prime number theorem for 
arithmetic progressions leads to the approximation 

3.1 h[^q)ip{x; q, l)\ ~ h(x) + . 

logx 

As a result 

S'(x,q) = B(x;q,R) - B(x;q,N) 
= \i[4>{q)^{x; q, R)] - li[(p(q)^(x; q, N)] + <f>(q)6(x, q) 
~<t>(q)[n'(x;q,N)-Tr'(x;q,R)}. 

The asymtotic equivalence in (13.11) holds up to the error term 2, p. 260] O(^f^), 
with 

R(x) = min (x e <> log 2 x, xe" ^ 71 ^) , a > 0, 

q = max K mod q (sup3?(jo), p a zero of L(s, «)). 
Let us now look at the statement GRH =>■ (i). Following [3l p 178-179], one has 

a) = —p- k{o)iIj(x,k) 

K mod q 
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Figure 3. The normalized regularized bias S'(x,13)/^/x versus 
the Chebyshev's bias <5(x, 13) at the prime champions of 5(x 1 13) 
(when <5(x, 13) > 0), and the curve at the prime champions of 
— 5{x, 13) (when 5{x, 13) < 0). The extremal prime champions in 
the plot are x = 263881 (with 5 = 123) and x = 905761 (with 5 = 
—40). The small dot curve corresponds to the fit of S'(x, 13) j yjx 
by 2/logx calculated in Sec. 3. 
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Figure 4. The normalized regularized bias 8'(x, ll)/y/x versus 
the Chebyshev's bias <5(x, 11) at the prime champions of S(x, 11) 
(when 6(x, 11) > 0), and the curve at the prime champions of 
— 5(x, 11) (when 5{x, 11) < 0). The extremal prime champions in 
the plot are x = 638567 (with 5 = 158) and x = 1867321 (with 
5 = —32). The small dot curve corresponds to the (very good) fit 
of 5'(x, 11)/ ' \fx by 2/ logs calculated in Sec. 3. 
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and under GRH 



ir(x) c(q,a) Jx i 
n(x;q,a) = — --—— + — — — ^ «(a^(x,«)+0( 



(f){q) (f>(q) logx <£(q)logx ' log 2 a; 

where Ko is the principal character modulo q and 

c(q, a) = -1 + #{1 <b<q:b 2 = a mod <?} 
for coprimes integers a and g. Note that for an odd prime q — p, one has c(p, a) — 

(?)■ 

Thus, under GRH 

TT(x;q,N)-7r(x;q,R) = ^^[^{c{q, R) - c(q, N)) 

(3.2) +e« mod M N )-~< R ))^ K )+°((^y 

The sum could be taken over all characters because Ko(iV) = kq(R). In addition, 
we have 

(3.3) i>(x;q,N)-iP(x;q,R) = -±-r £ [k(N) - k(R)]1>(x,k). 

k mod q 



Using (|3.2|) and (|3.3|) the regularized bias reads 

§'{x,q)r,- K '[x;q,N)-'K l {x;q,R) 

(3-4) = -cfo.AOl+O^ 

For the modulus q = 4, we have c(q, 1) = — 1 + 2 = 1 and c(q,3) — —1 so that 
<5'(x,4) = The same result is obtained for a prime modulus q — p since 

c(p, AT) = -1 and c(p, i?) = c(p, 1) = (±) = 1. 

This finalizes the proof that under GRH, one has the inequality ir'(x;q,N) > 
w'(x;q,R). 

If GRH does not hold, then using [2] lemma 2], one has 

B(x;q,a) — f2±(ar) for any £ < 9 q . 

Applying this assymptotic result to the residue classes a = R and a = N, there 
exist infinitely many values x = x\ and x = X2 satisfying 

B(xi;q, R) < —x{ and B(x 2 ;q, N) > x% for any £ < q , 

so that one obtains 

(3.5) B(xl;q, R) - B{x 2 ;q, N) < -x{ -x\<Q. 
Selecting a pair (xi,x 2 ) either 

B(xi;q,R) > B(x 2 ;q,R) 
so that B(x 2 ; q, R) — B(x 2 ; q, N) < and (i) is violated at x 2 , or 

(3.6) B(x 1 ;q,R)<B(x 2 ;q,R). 

In the last case, either B{x\\ q, N) > B(x 2 ; q, N), so that B(x\; q, R)—B(xi; q, N) < 
and the inequality (i) is violated at xi, or simultaneously 

B(xi;q,N) < B(x 2 ;q, N) and B(x 1 ;q,R) < B(x 2 ;q,R), 

which implies (|3.5|) and the violation of (i) at x = xi = x 2 . 
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To finalize the proof of l3.1[ and simultaneously that of ll.l[ one makes use of the 
asymptotic equivalence of (i) and (ii), that is if GRH is true => (ii) => (i), and if 
GRH is wrong, (i) may be violated and (ii) as well. 

Then, proposition ^ . 1 l also follows as a straigthforward consequence of proposition 

o 

□ 

11 

4. Summary 

We have found that the asymmetry in the prime counting function ir(x;q, a) 
between the quadratic residues a = R and the non-quadradic residues a = N for the 
modulus q can be encoded in the function B(x; q, a) [defined in introduced 
by Robin the context of GRH [2], or into the regularized prime counting function 
ir'(x;q,a), as in Proposition [3T] The bias in tt' reflects the bias in tt conditionaly 
under GRH for the modulus q. Our conjecture has been initiated by detailed 
computer calculations presented in Sec. [2] and proved in Sec. El Further work 
could follow the work about the connection of tt, and thus of tt' , to the sum of 
squares function r 2 (n) [ID] . 
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